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For a system of type-I neurons bidirectionally coupled through a nonhnear feedback mech- 
anism, we discuss the issue of noise-induced complete synchronization (CS). For the inputs to 
the neurons, we point out that the rate of change of instantaneous frequency with the instan- 
taneous phase of the stochastic inputs to each neuron matches exactly with that for the other 
in the event of CS of their outputs. Our observation can be exploited in practical situations 
to produce completely synchronized outputs in artificial devices. For excitatory-excitatory 
synaptic coupling, a functional dependence for the synchronization error on coupling and 
noise strengths is obtained. Finally we report an observation of noise-induced CS between 
non-identical neurons coupled bidirectionally through random non-zero couplings in an all- 
to-all way in a large neuronal ensemble. 
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I. INTRODUCTION 



Synchronous phenomena abound in nature and in our daily hves and have been studied from 
centuries past, right from Huygen's observations of synchronizing clocks. Various kinds of syn- 



chronous phenomena occur and have been identified (see for example 1 



13l|) (other references on 



the subject may also be found in Among these complete synchronization (CS) is one of 

the most interesting since the phase, frequency and amplitude of a subsystem all coincide with 
those of the other subsystem it is coupled to. It is seen therefore that in CS the trajectories of the 
coupled elements match exactly. CS is known to occur in identical systems, and was first demon- 



strated in chaotic systems in 



l|. Neurons and neuronal networks have been a subject of frequent 



theoretical and experimental studies 15]. Synchronization of neural activity has elicited a great 
deal of interest since it is believed that such phenomena enable cognitive tasks such as feature 
extraction and recognition to be performed. Hodgkin and Huxley classified neuron excitability 
mechanisms broadly into two classes 



Iq]: in type-II neurons the transition from a quiescent state 
to a periodically spiking state occurs through a Hopf bifurcation with a finite nonzero oscillation 
frequency. In type-I neurons, oscillations emerge through a saddle-node bifurcation on an invari- 
ant circle. As the bifurcation parameter changes, the stable and the unstable fixed points coalesce 
and then disappear, leaving a large amplitude stable periodic orbit. This is a global bifurcation 
and the frequency of the global loop can be arbitrarily small. Since axonal excitability patterns 
of mammalian neurons fall under the type-I class, it is but natural that this class has received 
special attention in the literature. Various observations have been made on type-I neurons, some 
prominent points of which are as follows. Equations for type-I neuronal dynamics can be reduced 
to the canonical normal form for a saddle- node bifurcation p^]. Repetitive firing occurs in the 
parameter regime when the system is in the close proximity of a saddle-node bifurcation on an 
invariant circle. Hansel et al [5| and Ermentrout have shown [61] that such neurons coupled via 



a certain class of time-dependent synaptic conductances are difficult to synchronize. Borgers and 
Kopell made further investigations of such coupled systems. In particular they discussed the effects 
of random connectivity on synchronization and the PING mechanism in networks of excitatory (E) 
and inhibitory (I) neurons both in the presence and in the absence of external noise 71, 8|]. 

In this work, we present some computer studies of generic type-I neurons coupled via synap- 



tic conductances such as those considered in 



7 



9[, which are governed by ordinary differential 



equations and which depend upon the outputs of the presynaptic neurons, and are subject to 
weak additive Gaussian white noise. We consider both excitatory-excitatory (EE) and inhibitory- 
excitatory (IE) bidirectional couplings and show that in certain regimes of the coupling constants 
and inputs, the system of coupled neurons shows complete synchronization (CS). The issue of CS in 
type-I neurons was not discussed in 7, 8]. As discussed in Section-2, largest Lyapunov exponents 
are shown to not adequately give information about CS in the system. We make an observation on 
the inputs to the neurons (and which are also modulated by the feedback in the system): we point 
out that the variation of the instantaneous frequency of the input received at each neuron with 
the instantaneous phase of the input it receives exactly coincides with that of the other neuron, 
in the event of complete synchronization of their outputs. It will be noted that in the presence of 
noise and feedback, this is not a trivial statement. We discuss the utility of this result in practical 
situations. 

In general, for EE synapses, our results indicate that when a common, externally applied con- 
stant input is used to perturb two bidirectionally coupled type-I neurons having identical coupling 
strength magnitudes and synaptic rise & decay times, weak noise induces them to exhibit CS upto a 
critical value of the coupling strength gc- For coupling strengths larger than gc, we find the system 
de-synchronizes through a power-law before locking on to a partially synchronized state for larger 
coupling strengths. We obtain a functional dependence for the synchronization error for neuronal 
outputs on coupling and noise strengths, in the regimes leading to partial synchronization. Such 



functional dependencies have not been reported in the hterature yet, to the best of our knowledge. 
In the noiseless case for identical EE neurons separated by different initial conditions, we observe 
that the antiphase states are stable in agreement with [6] and become completely in-phase in the 
presence of noise. For just two neurons with IE coupling, noise does not induce complete synchrony. 
In an ensemble of 200 non-identical neurons however, we find unexpectedly that noise-induced CS 
is possible with all-to-all bidirectional IE random couplings. 

II. COUPLED TYPE-I NEURONS IN THE PRESENCE OF GAUSSIAN WHITE NOISE 

In a system of n neurons, the activity of the ith neuron is described by a variable Xi which can 
be related to the membrane conductance. Its dynamics is represented by 

n 

Xi = qxf + li, /j = /3i + ^ Oij9jiSji , (1) 

where q denotes an inverse time constant for the membrane potential, Ij denotes its total input 
comprising of a constant external input /3j and the contributions from the presynaptic neurons, 
with i = 1, . . . n. Si j is the synaptic gating variable and represents the fraction of ion channels 
open in the jth presynaptic neuron, gij is the measure of the strength of the synapse from neuron 
i to neuron j; we have taken ga = 0. When Ij > this equation has no fixed points. Any 
initial condition tends to infinity in a finite time. To avoid this blow-up of solutions, a nonlinear 
transformation to new variables 6i may be made [3]: Xi = tan0j/2 which maps the real line onto 
a circle. Eqn.(2) then becomes 

^ = g(l - cos Oi) + (ft + ^ ajgjiSji{9j)){l + cos 6i) (2) 

i=i 

The point x = oo then gets mapped to the point 9 = ir and is interpreted as firing of a spike. We 
set = 1 and work in the parameter regime in which I/q « 1 so that the width of the spikes 
turn out to be in milliseconds as in real neurons 7|, |3] in these units. The Sij{t) evolve in time 



according to the differential equation which was considered in 7h9| 

dt Tij TR 

where Tij denotes the synaptic decay time and tr the synaptic rise time. The values of Sij always 
lie in the range to 1, reaching the maximal value when the neuron spikes. The synapse is an 
excitatory one if aj = +1, and —1 if it models an inhibitory synapse. Lest there be any confusion, 
we would like to clarify at the outset that when we refer to identical neurons, we mean neurons 
that have the same nature of synaptic coupling (EE) and have same coupling strengths 512 = g'21, 
have the same value of tr, tu = T21, and receive the same constant input /3i = /32, while their 
initial conditions differ by very little. 

Transmembrane voltage and neuronal firing can be affected by various sources of neuronal noise. 



but predominantly by synaptic noise la]. The synaptic noise itself occurs due to several factors, 
but chief among them is the synaptic bombardment at the inputs through the large number of 
neuronal connections, with each input spike adding a random contribution. We model this through 
an additive Gaussian white noise added to the neuronal input Ij. We study the dynamics of a system 
of two such neurons coupled bidirectionally as depicted in Fig.(l) and subject to Gaussian white 
noise ^(t) with the following properties: (^(t)) = , {(,{t)S,{t')) = 2a6{t — t'), where the stochastic 
variables are taken to obey Stratonovich calculus. Addition of Gaussian white noise ^ to eqn.(l), 
manifests as multiplicative noise in eqn.(2) because of the change of variables to 9, so that the 
equations now take the form 

d9 " 

^ = (1 - cos Oi) + il3i + Y, (^j9jiSji{9j) + + cos 9i) (4) 

Eqns.(3) for the Sij{9i) define the feedback regulating the activity of the jth neuron since Sij 
depends upon 9i which in turn depends upon 9j {i / j) via Sji{9j). The feedback increments 
or decrements the constant external input /3j received by neuron i. Thus the control parameter /j 
acquires time dependence through the dynamical variables. 



As in [lil, we consider the neuronal output to be described by the variable Ui = {1 — cos0j)/2 
as its time evolution pattern resembles that of a membrane potential in real neurons. This trans- 
formation maps the resting point Xi = corresponding to 9i = to Ui = 0, and the spiking point 
9i = IT to Ui = 1 via the relation Ui = xf/(l + xf). We choose to work with these variables as we 
get some new and interesting insights upon the dynamics underlying the phenomenon of complete 
synchronization. In terms of these variables, eqns.(4) and (3) become: 

N 



2{ui + {f3i + '^ajgjiSji + ^)(1 - Ui))y/ui{l 
^+exp(-2r?(l-n,))-^-^"'''^'^ 



Numerous studies in the 
various coupled systems 



iterature 
1 



(5) 



14 



rave reported the phenomenon of complete synchronization in 



2(>-l25l|. An adequately satisfying explanation of why and under 
what conditions CS can occur for systems with more complicated couplings, such as, for instance, 
that described by eqn.(5), is however, still lacking in our opinion. 



In particular, in 



21 



22l | the authors study noise induced CS in systems subjected to a common 
additive white noise and show that a necessary condition for CS is the existence of a significant 
contraction region in phase space. The systems studied in 



21 



22l | were Lorenz and Rossler systems 



which are far more amenable to analytical treatment than the equations above in eqn.(5). 
Since excitability in type-I neurons results from a saddle node on an invariant circle bifurcation, 
complete synchronization of two uncoupled neurons by common noise alone can be expected because 
of the existence of a contraction region close to the stable manifold of the saddle. On the other hand, 
when two such neurons are coupled together as in eqn.(5), the existence and nature of a contraction 
region would depend upon the eigenvalues of the Jacobian at the fixed points of the coupled system. 
However for eqns.(5), an analytical study becomes difficult since the Jacobian becomes singular 
at the fixed points. We therefore perform some computer studies on the system to learn more 
about the underlying dynamics. Moreover, since as we show below, Lyapunov exponents need not 



adequately give information about CS, we seek other explanations for occurrence of CS. 



As in 



21 



22l |. we define CS between the activities of neurons 1 and 2 as a vanishing value for the 
quantity {\ui — U2\) which is the synchronization error averaged over all iterations. 

Largest Lyapunov exponents for the system in eqn.(5) in the presence of noise for both EE 



and IE couplings were calculated following 



26H28t| and are shown in fig. (2). To incorporate noise 



29l43ll] was used. We note 



in the numerical calculations, the stochastic Runge-Kutta-4 method 
that in both EE and IE cases, the largest Lyapunov exponent Ai becomes more negative on the 
addition of noise. We observe also that for some intermediate values of the coupling constants 
(such as gi2 = 921 = 0.5), the Ai values could be larger, i.e., less negative, than those for lower 
coupling strengths (e.g., gu = g2i = 0.3). For EE coupling, Ai is almost always less than or equal 
to zero. In the case of IE coupling however, we find that for smaller noise-strengths and smaller 
couplings, for small (3, Ai fluctuates between positive and negative values in the presence of noise. 
This happens because of the oscillation of the bifurcation parameter (total input) between two 
regimes, depending upon the relative strengths of /3j and ajgjiSji, since aj = —1 for neuron j. 
Hence calculation of Lyapunov exponents may not adequately give information about CS or large 
windows of zero synchronization error, such as, for example, for the situation shown in Fig. (3), 
though they may certainly show the emergence of a definite order in the presence of noise and 
possible synchronization between the coupled units. 

Indeed, CS is expected to occur between identical systems and finding CS between non-identical 
oscillators would be unusual. As we describe later, however, we do find noise-induced CS in a 
system of 200 non-identical oscillators with random all-to-all couplings. Equally intriguingly, we 
find that there are parameter regimes where identical oscillators do not show CS at all, though the 
Lyapunov exponents shown in Fig. (2a) remain negative. Negative transverse Lyapunov exponents 
are widely accepted as characterizing CS, but their calculation for the system of 200 neurons, with 
a nontrivial feedback mechanism for phase resetting, as in eqn.(5) is a difficult task and we have 



8 

not attempted it here. In our case, even for n = 2, if we were to define new variables: z = 

and w = ("1+"^) ^ expressing the synchronization error dynamics through z in terms of z and w 

alone is not at all straightforward for the equations in (5). 

Moreover when each of the subsystems has a saddle-node-on-an-invariant circle bifurcation in the 
uncoupled limit, one could expect windows of intermittent firing patterns. In the coupled system, 
this is indeed observed in some parameter regimes, interspersed with large windows showing zero 
synchronization error, even in identical oscillators (EE case) for lower noise strengths (Fig. (3)). 
Hence we believe that negative largest Lyapunov exponents alone may not constitute conclusive 
proof for predicting noise-induced CS in coupled neurons. 

We therefore looked for other indicators which could help in understanding the mechanism of 
CS better in systems with feedback, such as in eqns.(5). We found one such simple indicator for CS 
in the context of the model under study and which we will now describe. The same methods and 
analysis should also hold for getting information on CS in any other system. Since synchronous 
activity is brought about by a common input or through mutual interactions and since these include 
components which are highly random, we study the instantaneous values of the sum total of the 
inputs received by each unit of a coupled system. We first set up a framework for this purpose and 
then provide a physical motivation and explanation for understanding CS through this indicator. 

III. INSTANTANEOUS PHASE - INSTANTANEOUS FREQUENCY VARIATIONS OF 

THE NEURONAL INPUTS & CS 



We construct the analytical signal [IJ, l2o|: Ti{t) = Ii{t) + iH{Ii{t)) = Si(t)e*^»W for the 
inputs h and similarly w{t) = u{t) + iH{u{t)) = i2(t)e'<^W for u{t) using Hilbert transforms. The 



instantaneous amplitudes and phases evolve according to: 



Ri{t) = Ri{l + Rf - 2Ri cos (f)i){Ri cos(Ci - A) - Ri^i cos(Ci -A + Pi) + Bi cosiipi + pi) + ^r} 
Mt) = 2^l^±^^^^^^{Ri sin(Ci - Vi) - RiBi sin(a - i^i + Pi) + Bi sin(Vi + Pi) + 

where Q^t) = arctan(-^^|^^^^), and ipi{t) = {Q{t) — (j)i{t))/2 denotes the instantaneous difference 
between the phase of the output of the ith neuron and that of the part it sends as feedback to the 
presynaptic conductance of the jth neuron — this feedback is a stochastic component of the input 
for neuron j. We have constructed the noise terms ^R{t) and C<^(t) from the analytical noise signal 
^{t) + i,yit) = m + iH{m)- 

(,R^ = {1 - Ri){^ COS 'ilJi + v sin ipi), = {1 - Ri){^ sin ipi + i/ cos ipi) (7) 

^Ri{t) and ^((,i{t) are periodically modulated by ipi which evolves according to the differential 
equation: 

ipi = I ={(-^1 - sin(Ci - Vi) - BiRi sin(Ci - + Pi) + Bi sin(V'i + Pi) 

'Ri{l + Rf-2Ri cos (f)i) 



+ - Ri sin 4ii[RiCOs{Q - in) - BiRiCOs{(i - ii + pi) + BiCOs{ipi + pi) + CrJ} 

= — ^ {{R^ -1)^/R^- sin (l)R} (8) 

2^{1 + Rf -2RiCos^i) 

The instantaneous phase therefore has no deterministic time scales and its drift and diffusion in 
time in the presence of noise is influenced by Bi and pi, the instantaneous values of the amplitude 
and phase respectively of the input, and also by the instantaneous amplitude Ri of the neuronal 
output. Bi{t) and pi{t) evolve as follows: 



1 g-2j7(l-i?j cos(/)j) . p. 

Bi{t) = —Bii 1 cos(2r7i2j sin^j) j H ^cosp. 



Tji TR ' Tji 

+ -(A + ^) cos((2r?i?,- sin<^,) - pi) e-2''(i-«. 



.J cos <t>j) 



TR ^ TR 

3 



g-2?7(l-_Rj cos(/>j) 

Pi{t) = sin(2r7i?jsin(?!)j) - — ^sinpj 

TR BiTji 

+ -^iPi + y ^) sin((2ryi2,- sin <Pj) - pi) e-^^^^-""^ ^°^'^^) (9) 

TRBi ^ TR 
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These show the effect of feedback on neuronal response. It is seen that CS between the outputs 
of neurons 1 and 2 occurs when the changes in the instantaneous phases and amphtudes of the 
two neuronal inputs exactly match each other. In other words, CS in the inputs to the neurons 
is required for the outputs to synchronize in phase, amplitude and frequency. This observation is 
in general not quite obvious since the system is nonlinear, and has a feedback mechanism which 
depends upon the outputs of the other neurons. 

In Fig. (4) we present the instantaneous-phase versus instantaneous- frequency plots of the inputs 
received by the two neurons. In all the numerous cases we studied for the coupled system for 
n = 2, we found that the signature of CS is the almost identical nature of these plots for the two 
systems that are in synchrony, be it with or without noise. On the other hand, the absence of CS 
gets reflected in the non-identical variation between instantaneous phase and instantaneous rates 
of phase-change of inputs to the two neurons, in the plots shown in Fig. (5). This is again true for 
both the noiseless as well as noisy cases. 
The system of equations (eqn.(5)) is of the form 

Ui = fD{ui,Si,Uj,t,l3 = O^ajQji = 0) + {Pi + ajgji,Si{uj,t))fn{ui) + Un{ui), (10) 

where 

fDiui,Si,l3^ = 0,ajgji = 0) = 2uf\l - Ui)^/\ Uui) = 2u]'\l - Uif/\ (11) 

and since feedback to each neuron through the synaptic coupling is oscillatory, in a Fokker-Planck 
description of the stochastic process, the probability distribution P{ui,U2,t) of the neuronal en- 



semble will not be stationary in the t — t- oo limit 



32( 1 . In certain regimes of the noise strengths and 



coupling constants, where competing contributions from the drift and diffusion terms would make 
the noise-averaged difference in outputs {{ui — U2)) zero, CS occurs. 

Physically, CS is brought about through the following sequence of events. Addition of a small 
amount of noise increases the decay time of the synaptic conductances Sij gradually, and eventually, 



11 

lowers their minimum to zero. This delays the onset of the next Sij peaks, and hence the input li 
to each neuron at any further instant of time. Increasing the noise strength further increases the 
decay time of Sij. The periodically maximal values of the inputs thus take longer to arrive at the 
neurons and this becomes visible in the neuronal firing pattern as departures from the previous 
(noiseless) values of the phase differences between the neurons, and those of their output differences 
(■ui — U2)- The instantaneous values of the phases and the rate at which they change in time, i.e., 
the instantaneous frequencies of each neuron, is determined by the strengths aijgij of the synaptic 
inputs it receives and the noise strength for any given set of /3,TR,Tij,r]. Hence for given initial 
conditions and different amplitudes for ui &l U2, it would be reasonable to expect CS to occur when 
the following is satisfied for the inputs to the neurons: the variation of instantaneous values of the 
frequencies with instantaneous phases of the input for neuron 1 matches with that for neuron 2. 
This results in the instantaneous values cpi and <pi of neuron 1 changing in step with (p2 and (p2 
respectively of neuron 2. This forces the amplitudes of neurons 1 & 2 to become identical with each 
other, since otherwise both conditions (pi — (p2 and (pi = (p2 cannot be simultaneously maintained. 
Thus CS results. 

A striking feature of all these plots is their strange, flame-like structure. The flame shape is 
reminiscent of canards that are typically associated with systems exhibiting relaxation oscillations. 
Indeed, from the common factor occurring in the inverse square root on the right hand side of 
eqn.(8), it is apparent that ipit) does evolve on a time scale different from that for R{t), (p{t),B{t) 
h p{t) in eqns.(6) &; (9), though it is a dynamically varying time scale, determined also by the 
stochasticity of the system. The separation of time scales is in fact manifested in the time series 
for the neuronal inputs which show relaxation oscillations (Figs. (4) &; (5)). This gives the input 
instantaneous phase-frequency curve its characteristic shape whenever noise is introduced. This 
is indicative of some order emerging in the phase. In fact, noise-induced phase synchronization 
can be demonstrated even in the IE system by interpreting the instantaneous phase differences 



12 



$ = — 02 in a statistical sense as in 



2ll |. The distribution -P($), of cyclic instantaneous phase 
differences, ^ £ [— vr, vr] is shown in Fig. (6). Preferred phase differences between the 2 neurons 
manifest as peaks in <I> which become sharper and taller with increasing noise-strengths - a clear 
indication of noise-induced phase coherence. 

IV. NOISE-INDUCED CS & PARTIAL SYNCHRONIZATION IN EE &: IE SYSTEMS 

In Fig. (7a) we have plotted the synchronization error {\ui — U2\) for two coupled identical (EE) 
neurons and for coupled non-identical (IE) neurons, as a function of the coupling strength for 
different noise-strengths. We have considered here the special case of /3i = /32 = 0, and each 
neuron receives inputs only from the other neurons; i.e., we look at the effect of feedback. Results 
and analysis for non-zero /3j are presented elsewhere jssj]. We see that although CS is expected 
between identical oscillators, it does not happen for the noiseless (deterministic) case. Increasing 
noise-strength brings down {\ui — U2\) and indeed takes it down to zero for certain ranges of the 
coupling constant. 

For the identical EE coupled neurons, when feedback constitutes the only input (there being no 
other explicit input), as in Figs. (7a, 7b), then beyond a maximal critical coupling constant strength 
Qc, the system gets de-synchronized, with {\ui — U2\) 7^ 0. For g > gc, there then exists a regime 
in the de-synchronized system where noise strength a still plays a role in determining the output. 
In this transition regime, for g > g^ we find that the difference in outputs — U2\)trans depends 
upon g through the expression 

{\ui - U2\)trans = {O - 9cY^^ " Z^cr, (12) 

where Va is a constant that depends upon the noise strength a. 

The critical coupling strength g^ depends upon the noise strength as well and we find that it varies 



13 



as 

£/c ~ oaA, (13) 

where a = 1.1. At very high g the system gets locked to a partiaUy synchronized state, with 
(|iti — «2|) approaching a constant value, wherein noise strength no longer influences the difl'erence 
in the outputs of the neurons. In Fig.(7a), this constant value approaches 0.5. For the curves in 
Fig.(7a,b), for the entire regime following the beginning of desynchronization, we find a functional 
dependence on the coupling constant g given by an equation of the form 

{\ui-U2\)a ^ a{cF)-h{(T)g- (for 5 > 5c), (14) 

where a{a) and 6(17) depend on the noise-strength through 

0.59 0.0181 

so that 

0.59 0.0181 1.4641 om 

dm - «2 Da ~ ^-^0:53-5-^3—^ ■ (for 5>5e)- (16) 

This expression is plotted on the numerical data points in Fig. (7b). A rigorous theoretical treatment 
of the system needs to be done in future studies to establish these relations for the synchronization 
and desynchronization transitions through a Fokker-Planck approach. 

In the IE case (Fig. (7c)), we see that noise-induced CS does not occur for two coupled neurons; 
however, there is partial synchronization (by this we mean that {\u\ — 1*2!) = constant 7^ 0) since 
(1^1 — U2\) gets locked to a finite, non-zero value ('^ 0.6) for large g. Further, increasing noise- 
strength increases rather than decreases {\ui — it2|), in the region before partial synchronization, 
in contrast to the observation for the EE case. We were unable to achieve CS in the IE case for 
two coupled neurons, even in the presence of noise. 

However, for an ensemble of 200 coupled theta neurons of which 150 neurons are excitatory 
and 50 inhibitory, we obtain very different results. In this simulation shown in Fig. (8), each 
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neuron receives the same input /3j = 0.1 and there is all to all random coupling with the coupling 
strengths varying between 0.05 and 0.1 and with different initial conditions. The excitatory neurons 
are shown in red while the inhibitory ones are in green. On introducing Gaussian white noise 
of strength a = 5.0 into the system we observe synchronous phenomena emerging between the 
excitatory and inhibitory neurons. The interesting thing to note here is that not only do most 
of the inhibitory neurons fire in synchrony with other inhibitory neurons, but that also most of 
them are synchronized with the excitatory neurons. This kind of noise-induced near-CS in coupled 
type-I neurons with random non-zero coupling strengths has not, to our knowledge, been reported 
previously in the literature. A detailed study of this situation is under way to explain the observed 
synchrony for n = 200 and will be reported elsewhere: it is beyond the scope of the present work. 
Interestingly, in the literature, spatiotemporal synchronization has been shown to occur in networks 
of coupled chaotic maps with varying degrees of randomness in the coupling connections I4] . 



V. CONCLUSION 

We have studied the issue of noise-induced CS in coupled type-I neurons, a class of neurons 
that are especially important since mammalian neurons fall under this category. We find that 
complete synchronization between any two neurons is signalled whenever the variations of the 
instantaneous input phases (p versus the instantaneous frequencies (p of the inputs of the neurons 
being studied are identical. We point out that such identical (p-cj) plots of the neuronal inputs would 
be a signature of CS between the neurons. This suggests the possibility of producing completely 
synchronized outputs of coupled systems having feedback mechanisms in the presence of noise, by 
ensuring that the plots of the instantaneous frequency versus instantaneous phase of the inputs to 
the subsystems are identical. That this is a significant point will be appreciated when one recalls 
that this condition is required in the continued presence of noise and feedback. When monitoring 
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or control of neurons is required to be done in a living organism in the case where synchronized 
neuronal output is required at another, inaccessible spot through a given external input, this result 
becomes important. 

Though CS is expected between coupled, identical neurons, we find that CS occurs only upto 
a critical value gc of the coupling constant g for a given noise strength beyond which the system 
desynchronizes again, and then for large g gets locked to a partially-synchronized state. We find 
that the critical coupling strength gc depends upon the noise-strength through a power law. For 
g greater than the critical value gc, from the transition regime through upto the onset of partial 
synchronization, we find a functional dependence of the noise-averaged output difference on g and 
the noise strength a, given by eqn.(14). For a larger ensemble of 200 neurons, we find unexpectedly 
that non-identical neurons can show near-complete synchronization. 

Since type-I neurons model axonal excitability patterns in mammals, the results presented 
here would be useful in the study of synchronous mechanisms underlying the neural code. As 



an immediate application, we believe our results would be usefu 
observations reported on cat & awake monkey visual cortex 34- 



in explaining the experimental 



38l | which show synchronization 



of neuronal activity with a single stimulus, and which disappears when activated by different, 
independent stimuli. Since a single stimulus would correspond to neuronal inputs of identical 
amplitudes, instantaneous phase and frequency, this is actually the same scenario that we have 
found for CS to occur. It is also likewise clear why CS was experimentally observed to vanish 
on the activation of different stimuli, since the necessary conditions of identical input amplitude, 
phase and frequency are no longer met. 
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Figure 1. Bidirectional coupling between two theta neurons. Left: E-E coupling, Right: I-E 
coupling. 
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Figure 2. Lyapunov exponent in the presence of noise for: (at top) EE coupling 
and (at bottom) IE coupling. In both cases, the parameters which we used are 
Pi=h = 0.1, ri2 = T2i = 2.0, TR = 0.1, rj = 5.0. 
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Figure 3. Windows of zero synchronization error interspersed with intermittent bursts of firing 
in coupled identical EE neurons: 1712 = 521 = 0.3, T12 = T21 = 2.0, tr = 0.1, /3i = /32 = 0.1, 
a = 0.2655. 
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Figure 4. Instantaneous phase-frequency "flame" plots of stochastic input Jj = A + ^^x ajgjiSji 
to neuron i, for EE coupling, for (a) cr = (top) and (b) a = 0.09 (bottom). In these plots, 
Pi = ^2 = 0.0, showing the effect of feedback. The extended panels below each set of flame plots 
show the corresponding time-series for neuronal input (at left) and the difference (ui — U2) in 
neuronal outputs (at right). Parameter values are: gi2 = g2i = 0.3, T12 = T21 = 2.0, tji = 0.1. 
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Figure 5. More instantaneous input phase- frequency "flame" plots with corresponding time-series 
of neuronal input, row 1: EE, a = 0.0; row 2: EE, a = 0.301; row 3: IE, a = 0.0; row 4: IE, 
a = 0.301. Parameters used are:5fi2 = 521 = 0.3, T12 = T21 = 2.0, /3i = /32 = 0.1, tr = 0.1 . 
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Figure 6. Distribution of instantaneous phase-differences for neuron pair with IE coupling, 
different noise strengths. Here, gi2 = 521 = 0.3, ri2 = T21 = 2.0, tr = 0.1, = 5.0, /3i = /32 
0.0, «! = +1,02 = —1. 
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Figure 7. Noise-induced synchronization in coupled Type-I neurons (a) EE case; transition 
from synchronized to partially synchronized state when g > Qc (b) EE case, as in (a); solid lines 
correspond to Eqn.(16). (c) IE case; noise-induced CS is absent but there is partial synchronization 
as the system gets locked to {\ui — U2\) ~ 0.6 for large g. For both (a) & (b): ri2 = T21 = 2.0, 
(3^= (32 = 0.0, TR = 0.1, initial conditions 9i = 0.0, 6*2 = 0.01, S12 = S21 = 0.0 
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Figure 8. Simulation of 200 bidirectionally-coupled theta neurons (150 EE (in red) and 50 IE (in 



green)) with random all-to-all couplings; top: a = 0, bottom: a = 5.0. 



